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Abstract 

We discuss the issue of duality violations in hadronic tau decay. After in- 
troducing a physically motivated ansatz for duality violations, we estimate their 
possible size by fitting this ansatz to the tau experimental data provided by 
the ALEPH collaboration. Our conclusion is that these data do not exclude 
significant duality violations in tau decay. This may imply an additional sys- 
tematic error in the value of a s (m T ), extracted from tau decay, as large as 
5a s {m T ) ~ 0.003 -0.010 . 



1 Introduction 



The hadronic decay of the tau lepton is a particularly suitable process for studying QCD 
interactions, because its mass is heavy enough to produce hadrons in the final state, 
without the complications which arise with hadrons in the initial state. However, since 
the mass of the tau is not much larger than a typical hadronic scale of order 1 GeV, it 
was initially unclear whether a perturbative QCD analysis of this process would not be 
ruined by nonperturbative effects beyond any systematic control. The pioneering work 
of Ref. pQ showed, however, that perturbation theory, augmented with the Operator 
Product Expansion (OPE) [2], can indeed be the right tool for understanding strong 
interaction aspects of hadronic tau decay. Among these, the most important ones are 
perhaps the determination of the coupling constant a s , the light quark masses, and the 
vacuum condensates [3]- 

In fact, the determination of a s from tau decay is now one of the most accurate ones. 
The level of accuracy is so high that subtle effects, which were completely negligible 
before, may now start to play a significant role in the precision of the result. Since 
the value of a s is of fundamental importance for our understanding of QCD and the 
Standard Model, it is obvious that a good control of all systematic effects is, now more 
than ever, mandatory. 

References [3 El [71 [8] find the following values for a s after analyzing the same 
experimental data obtained by ALEPH [I]: 



Note that, if the errors are taken at face value, these four results are not fully consistent 
with one another. The main difference between the results of Refs. [5] and [8] is due 
to the treatment of perturbation theory (the use of so-called Contour-Improved (CI) 
vs Fixed-Order (FO) resummation prescriptions, respectively [§]). This difference has 
been included in the systematic error in [Bj. The difference between the result of 
Ref. [7J and the other three is mainly due to the use of different weight functions in 
the finite-energy sum rules. Different choices for these functions distribute differently 
the relative weight between the perturbative and the OPE contributions, as extracted 
from the spectral function. Although in an exact treatment these choices should not 
matter, in practice they do since truncations of both the perturbative series and higher 
orders in the OPE have to be applied. Ref. [7J also analyzes the effect of some of these 
higher orders in the OPE. 

Another interesting "mismatch" is seen in the analysis of the ALEPH data carried 
out in Ref. [5], which quotes the following values for the gluon condensate: 




0.344 ± 0.005 exp ± 0.007 th [5J , 

0.332 ± 0.005 cxp ± 0.015th [fij , 

0.321 ± 0.005 cxp ± 0.012 th [7J , 

0.316 ± 0.003 cxp ± 0.005 th [H] . 



(1.1) 



a. 



s 



(GG) 



(-0.8 ± 0.4) x 10~ 2 GeV 4 , 



(1.2) 



7T 



Vector 



7T 



s 



(GG) 



Axial 



(-2.2 ± 0.4) x 10~ 2 GeV 4 , 
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where the subscript denotes whether the condensate has been determined using the 
vector or the axial-vector spectral function. Although these two quantities should be 
one and the same in QCD, it is clear that the two values in Eq. (11.21) are not compatible 
with each other. Of course, the discrepancies (11.11) and (11.21) would not be there if either 
a) one could argue that a specific analysis is to be preferred on theoretical grounds, 
or b) if the errors had been underestimated because of yet other systematic effects. 
Further determinations of the OPE condensates may be found in Refs. [10]. 

From the theoretical point of view, one could think of (at least) three sources of 
possible systematic errors, namely the truncation of the perturbative series in powers 
of a s , the truncation of the OPE, and the contribution from Duality Violations (DVs). 
In fact, at some deeper level, the three sources must be related. First, the purely 
perturbative contributions can be viewed as the contribution from the dimension-zero 
"condensate" (i.e., the unit operator) to the OPE. Second, the lack of convergence of 
the perturbative series in a s calls for the existence of some nonperturbative contribution 
which, through renormalons, is related to the OPE condensates. Finally, the lack of 
convergence of the OPE is, in turn, at the origin of DVs. 

It is far from clear that the different values for a s quoted in Eqs. (II .ip should be 
due to DVs, since the first two error sources mentioned above can potentially explain 
this difference by themselves. However, in this work we will concentrate on DVs as a 
further possible source of error, because it has been much less explored than the other 
two. 

In QCD analyses of tau decays, the OPE plays a central role. However, although 
this expansion is of such fundamental importance, the properties of the OPE in QCD 
are not known. In particular, although it is suspected that the OPE is asymptotic, 
this is not known to be true, and, even if it is indeed asymptotic, we do not know 
whether it is summable or not, or how it behaves along different rays in the q 2 complex 
plane. For an asymptotic expansion, it is of course important to estimate its intrinsic 
theoretical error. Relegating a more precise definition of this error to the discussion in 
the next section, here we just mention that this intrinsic theoretical error is normally 
associated with DVs. 

In practice, one can perhaps get a feeling of the systematic error in the perturbative 
expansion by comparing different orders in the a s expansion. Likewise, the systematic 
error in the condensate contribution may also be assessed from the comparison between 
different orders in the OPE. However, we do not know of any systematic approach for 
studying DVs. Therefore, although the disagreement between, e.g., the results quoted 
in Eq. (11.21) may point towards the fact that "this method may approach its ultimate 
accuracy" in tau decay [5], there is currently no method to estimate this accuracy from 
first principles in QCD. Hence, there appears to be no other way to make progress than 
to resort to models. 

Based on tau data, the only attempt to date at estimating the error from DVs is 
the one made in Ref . [5] ; their claim is that DVs effects are completely negligible at the 
tau mass. This result was obtained from the V+A spectral function with the help of a 
physically motivated model, previously used for studying aspects of DVs [T ^ [T5 | [T6l [T7] . 
In the present work, we reanalyze this estimate studying separately the V and A 
correlators. As we will see, our conclusion is different: based on tau data alone, DVs 
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may be larger than those found in Ref. [5]. The inclusion of e + e~ data beyond the r 
mass, with some extra assumptions to be detailed below, may also be used to further 
constrain the size of DVs but, even in this case, their impact could be comparable to 
the presently quoted systematic uncertainties. Therefore, the possible effect of DVs 
should be taken into account. 

This work is organized as follows. In Section 2, we present the theoretical back- 
ground needed for the discussion of DVs in tau decay, while in Section 3 we introduce 
and motivate our ansatz for DVs. The actual fits to the experimental tau data are 
described in Section 4. Section 5 is devoted to an estimate of the impact of these 
fits on the determination of a s and, in Section 6, we assess the changes imposed on 
this estimate if one also considers recent e + e~ data beyond the r mass. Finally, in 
Section 7, we offer some conclusions and prospects for future work. 



2 Duality violations and r decay 

Let us start our discussion of hadronic tau decay by defining the correlators 

n^(g) = iJd*xe**(0\T{J^x)Jt(0)}\0) 

= W ~ q 2 9n ImII« (a) + q^U^ A (q 2 ) , (2.1) 

in which JY(x) = u{x)^^d{x) and J^(x) = u(x)^ ^ 5 d(x) . The ratio of the decay 
widths of the tau lepton to nonstrange (vector and axial-vector channel) hadrons and 
the decay width to electrons, denoted as 

s r[T-~(W,] 

r r- -► v T e-v e \ 



can be expressed as 



RY' A = 12nS EW \V ud \ 2 l ^li- 



ds 



ml V m 2 . 



s 

mi 



1 + 2—1 Imn^(s) + lmU^ A (s) 



(2.3) 

where the factor Sew — 1.0201(3) [llj accounts for small, known electroweak correc- 
tions, and V u d is the corresponding entry in the CKM matrix. 

The spectral functions Im IIy^(s)0 appearing in Eq. (12.31) have been measured most 
accurately by ALEPH [4J, and the results are included in Fig. [3] below. Since at present 
these spectral functions cannot be directly calculated from QCD, what is being done 
in practice, following Refs. [1] [T2J, [13] , is to take advantage of the analytic properties of 
the functions Uy\(q 2 ) in the complex q 2 plane, depicted in Fig. [TJ We write (dropping 
the superscript ^ from now on) 



dsP(s)-lmU v>A (s) = ~ <b dq 2 P(q 2 )U VtA (q 2 ) , (2.1) 



-'JmlTy^s) vanishes in the isospin limit, which we assume, and Imn^(s) is dominated by the 
pion pole. 
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Figure 1: Analytic structure ofHv,A(q 2 ) in the complex q 2 plane. The solid curve shows the 
contour used in Eq. \2.1$ . 



where P{q 2 ) is an arbitrary polynomial which can be chosen at one's convenience. In 
the application to tau decay one sets the radius so = m 2 and, under the assumption 
that this scale is larger than the scale at which the OPE sets in, it makes sense to 
rewrite the right hand side of Eq. f)2.4p as 



1 

2ni 



I dq 2 P{q 2 ) |lT^V) + A VA (q 2 )} , (2.5) 

J\q*\=8 1 J 



where, by definition, 



A VA (q 2 ) = U VA (q 2 ) - n° P A V) , (2.6) 

and then approximate Ay,A{q 2 ) ~ 0, thus neglecting the second term in Eq. (12.51) . To 
the best of our knowledge, in all analyses of tau decay to date, except for Ref. [5], this 
approximation has always been assumed. As we already mentioned in the introduction, 
Ref. [5] does discuss the possibility of DVs in tau decay, but only for the combined 
V+A correlator, and concludes that DVs are negligible and do not affect the systematic 
errors from other sources, cf. the first result quoted in Eq. (ll.ip Fl 

Here we will not assume this approximation. On the contrary, we will take non- 
vanishing contributions to Av,A(q 2 ) as defined in Eq. (12.61) as the definition of DVs. 
Consequently, if we define 



2>f>o) = -7T-f dq 2 P(q 2 ) A VA {q 2 ) , (2.7) 
2m J\ q 2 l=S0 

as a measure of these DVs in tau decay, Eq. (12.4ft becomes 

f° ds P{s) -Imn y , A ( S ) =~<f dq 2 P{q 2 ) Il° P A E (q 2 ) + V [ ^(s ) . (2.8) 
2 This conclusion will be reassessed in Section 5. 
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Re q 2 



Figure 2: Contour used in the derivation of Eq. )12.9\) . With the assumed exponential decay, 
Av,a(q 2 ) vanish on the dashed circle with infinite radius. 



Note that, if the OPE were a convergent expansion and Sq were within the radius of 
convergence, one would get Ay^(q 2 ) = and thus DVs would be absent. 

Although the precise convergence properties of the OPE in the complex plane are 
not known, we do know that this power expansion cannot be convergent but only 
asymptotic (at best). A convergent expansion around \q 2 \ = oo in the complex plane 
defines an analytic function on a complete annulus around the origin, which must, 
therefore, also include the Minkowski axis. However, this is contradicted by the exis- 
tence of the physical cut along the Minkowski axis, which shows that the expansion 
cannot be convergent for any sq- 

In this work we will assume that the OPE is an asymptotic expansion in q 2 and 
explore the possible consequences. As an asymptotic expansion, the OPE will have 
a region of validity (the so-called wedge of asymptoticity [18]) which will naturally 
exclude the Minkowski axis. This expectation is strongly supported by the large- 
N c limit, for which the physical cut becomes an infinite set of poles which are not 
reproduced by the OPE. We remark that, even if the function P(q 2 ) is chosen such 
that it vanishes on the Minkowski axis (as is the case for the polynomial appearing in 
the integrand of Eq. (12.31) ). that does not guarantee that the DV function T>y A (so) will 
exactly vanish. 

Accepting that the OPE is asymptotic, it is natural to assume a momentum de- 
pendence of the DV function Ay.A{q 2 ) with an exponential decay at large |g 2 |Jfl Then, 
using the contour depicted in Fig. [2] and the exponential decay of Ay^iq 2 ) on the 
dashed circle at infinity, one obtains a simpler expression for T>y A (so) [HH [T7] : 

/•OO 1 

V [ v ] A (s ) = - ds P(s) -ImA v>A (s) . (2.9) 

This relation is interesting for two reasons. It is useful, since it expresses the DV 
contribution to Eq. (jZEl . V [ y\ (so), as an integral along the Minkowski axis, in terms 

3 This is the analog of the e _1//a -type error one expects in an asymptotic expansion in powers of a. 
Analyticity requires, however, the exponent to be angle dependent; see, e.g.. Sec. Ill B in Ref. [16] . 
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of the spectral function, and it can thus in principle be determined from experimental 
data (as we will do below). However, this requires an extrapolation from values of s 
for which experimental data are available to infinity. This shows explicitly the inherent 
difficulty present in any quantitative evaluation of DVs [19]. 

3 An ansatz for Duality Violations 

To make progress one clearly needs information about the functions ImAy^(s). How- 
ever, as emphasized above, it is not known how to get this information from QCD. 
Therefore, there is no other option but to adopt an ansatz which is theoretically sen- 
sible and which is, of course, not ruled out by existing data. In other words, we look 
for a consistent picture of how duality violations might occur in QCD, and how they 
may affect the existing determinations of a s (and, consequently, also quark masses and 
condensates). This way we are able to make an educated guess about the possible size 
of the systematic errors due to DVs as allowed by existing data. 

We will parametrize the duality- violating part of the spectral functions as 

-lm/\ V)A (s) = k V)A e~ lv - AS sin {a VA + Pv,a s ) > (3- 1 ) 

71 

The exponential decay in this expression is inspired by our assumption that the OPE 
is an asymptotic expansion, and can be understood as representing the finite width of 
resonances [H] . The oscillatory behavior is what one expects in a spectral function with 
resonances exhibiting some kind of periodicity. The precise form chosen here is that 
obtained if the vector and axial-vector resonances lie on Regge daughter trajectories. 

Furthermore, a glance at Fig. [3] reveals that both the vector and the axial- vector 
spectral functions cross perturbation theory around s = 2 GeV 2 . A natural interpreta- 
tion is that the corresponding DV functions ImAy^g 2 ) have a zero near this energy. 
Such a zero is much harder to understand without DVs since, in this case, it would 
have to be due to a cancelation between different orders of the OPE, signalling a pos- 
sible breakdown of the OPE. Although we cannot exclude such a breakdown at this 
rather high scale for minkowskian momentum, much evidence is consistent with the 
OPE being valid to a significantly lower scale for euclidean momentum, whether this 
scale is associated with the gluon condensate or with the quark condensate. Such a 
hypothetical breakdown is also not seen in the purely perturbative contribution, which 
remains small and essentially flat all the way down to 1 GeV 2 . Furthermore, as we will 
see in the next section and the Appendix, the OPE yields very small contributions to 
the V and A spectral functions from the condensates of dimension four and six down 
to 1 GeV 2 . 

Since the width of a resonance is a 1/N C effect, one expects the exponent 7 to 
be relatively small. The reason is that it has to go to zero when N c goes to infinity, 
since the exponential suppression shown in Eq. (13.1 1) has to disappear in favor of an 
unsuppressed contribution representing the isolated, infinitely narrow peaks present in 
any spectral function in QCD in the strict N c —>■ 00 limit. For the other parameters we 
expect k to be of order N c , and a and f3 are expected to be of order one in this limit. 
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Figure 3: Different combinations of spectral functions (V , A, and V ± A) (data points in 
black) compared with fixed-order perturbation theory (green flat line) and the result of our fit 
to Eq. (EJPj ( re d curve). 
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Although our arguments are admittedly rather heuristic, it is reassuring that a 
mathematical model realizing all the above features actually exists. In this model, first 
proposed in Refs. [T4"l [To] , resonance widths are introduced in a manner compatible 
with the analytic properties of the correlators (12 .11) . with masses following Regge theory. 
This model is physically well motivated, and has been used for a variety of studies of 
DVs for both light and heavy quarks [151 EEl CEO E]- Our ansatz, Eq. A3. If) . actually 
corresponds to the asymptotic behavior of this mathematical model for large s. By 
considering the approximation (13.1 p with a priori unknown values for the parameters 
a,/3,7 and k, we expect to describe the generic features of DVs at large s, avoiding 
some of the more model-dependent details. We note, however, that the ansatz (13.11) 
is a very good approximation to the full mathematical model already at relatively low 
values of s [16J. 

4 Fit to tau data 

To find out whether the ansatz (13. ip is compatible with the tau experimental data, we 
have fitted the ALEPH spectral functions in both the vector and axial-vector channels, 
i.e., we fitted the spectral data with the functions 



The function p(s) above contains the purely perturbative corrections in powers of a s 
up to 0{al)% These perturbative corrections differ depending on whether one uses CI 
or FO perturbation theory and the expressions for both treatments have been taken 
from Ref. [8], where one can also find a discussion of their relevance for tau decay. 

In principle one should also include in Eq. (14. ip the contribution from the con- 
densates which, away from q 2 = 0, can only come from the logarithms in the Wilson 
coefficients. However, the contribution from the operators of dimension four and six 
turns out to be numerically very small for reasonable values of the quark and gluon 
condensates. Our estimates for these contributions are relegated to the Appendix. 
Condensates of dimension eight or higher will not be considered. 

As already mentioned in the previous section, the ansatz (13.11) is to be understood 
as an estimate of DVs at large values of s. Since, in practice, we do not know how large 
s would have to be, we have done our fits in a window s min < s < m 2 , varying the 
lower end, s min . In other words, in the phenomeno logical approach we will be taking, 
the precise meaning of large values of s will be set by the quality of the fit. The results 
of the fit are very insensitive to the value of s m , n , between a maximum value of s m i n 
of ~ 1.8 GeV 2 (above this value the experimental errors are too large for the fit to be 
meaningful), and a minimum value of s m i n of ~ 1.1 GeV 2 . Below this value there is a 
sharp increase in the chi-squared per degree of freedom, showing a deterioration of the 
fit, probably due to the tail of the p meson (in the vector channel). Consequently, we 
have chosen 1.1 GeV 2 < s < m 2 as our fitting window. 

4 The term of 0{a b s ) in not known exactly and we have used the estimate provided by Ref. [5]. 
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(4.1) 
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The result of the fit to the vector spectral function is then: 



k v = 0.018 ±0.004, 

7v = 0.15 ± 0.15 GeV 

a v = 2.2 ±0.3, 

Pv = 2.0 ±0.1 GeV" 2 

x 2 io 



dof 79 
and that to the axial-vector spectrum: 

k a = 0.20 ±0.06, 

7a = 1.7 ±0.2 GeV~ 2 . 

a A = -0.4 ±0.1 , 

P A = -3.0 ± 0.1 GeV -2 

X 2 17 



— ~ 0.13 , (4.2) 



dof 78 



0.22 . (4.3) 



Even though the perturbative term p(s) in Eq. ( 14. ip is different in the CI or FO 
resummation schemes, the results for the fit parameters are, within the quoted errors, 
insensitive to this difference, or to the initial value for a s , which we took to be a s (m T ) = 
0.3156. Our fits are also rather insensitive to an increase of the value of s m , n ; the main 
effect is an increase of the errors quoted in Eqs. (14.2|4.3[) . 

Since our lower end is at s — 1.1 GeV 2 , one may worry about the use of perturbation 
theory at such low scales. However, there is nothing in the perturbative expansion 
signalling a breakdown of the approximation, as the nearly horizontal line in Fig. [3] 
clearly shows. Also, one may object that perturbation theory and the OPE should 
not be used for describing a spectral function on the time- like axis, Eq. ( 14.11) . i.e., 
that "local duality" does not work. However, we want to emphasize that local duality 
does not work precisely because DVs are not included, and that DVs, parametrized 
for example as in Eq. (14.11) . amount to adding the necessary contributions which make 
up for the missing piece, therefore making local duality exact (c/. Eq. (12.81) ). Our 
analysis thus amounts to the assumption that Eq. (14.11) gives a reasonable account of 
DVs above s = s min . 

Looking at the values obtained from the fits (14. 2114. 3p . one immediately realizes 
that the oscillations in the axial-vector channel will be more damped than in the 
vector channel as a result of a much larger exponent, i.e., ^a 3> Jv- While we have no 
understanding of why this should be so in QCD, we simply observe that the difference 
in the values of these parameters reflects the difference in the corresponding spectra as 
measured by ALEPH. 

Fig. E] shows the comparison of the fits to the different combinations of experimental 
spectral functions. As one can see, all fits are very good, having a very small \ 2 /dof . 
Our fits are fully correlated fits, i.e., correlations among errors were taken into account 
in both fits using the pertinent covariance matrix provided by ALEPH. The small 
values of x 2 /dof that we find might suggest an underestimation of these correlations 
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(or an overestimate of the diagonal errors), although it is difficult to know whether 
this is actually the case. In fact, an uncorrelated fit using only the diagonal entries 
of the covariance matrices leads to results which are similar to those in Eq. (14. 2114. 3|) , 
in particular for the vector spectral function. We remark that our analysis is not the 
only one to get such small values for \ 2 /dof . The analysis of Ref. [5] finds a value of 
X 2 /dof = 0.07 in their OPE fit to moments such as Eq. (12.31) in the vector channel. 

As can be seen from Fig. [31 the vector and axial-vector spectral functions oscil- 
late around perturbation theory in the window 1 GeV 2 < s < m 2 T with almost com- 
plete anti-correlation, with a common "oscillation node" around 2 GeV 2 . Therefore the 
combination V + A agrees with perturbation theory better than V or A individually. 
Furthermore, the fits of Refs. [H [5j suggest that there is a partial cancelation in the 

V + A combination of the contribution to the contour integral in Eq. H 2 . 5 1) from both 
the dimension-six and dimension-eight condensates in the OPE. These features have 
led people to believe that the combination V + A is a better correlation function than 

V or A, because it is somehow protected from all kinds of nonperturbative effects, 
including DVs. However, as the third panel on Fig. [3] shows, the partial agreement 
of perturbation theory (green flat line) with the experimental V + A data does not 
guarantee the absence of DVs. In fact, the results of our fits reported in Eqs. (14.2114.31) 
reproduce the experimental data better than perturbation theory and, consequently, 
constitute a clear indication that nonzero DVs may exist in this spectral function as 
well. Furthermore, because DVs are suppressed in the axial-vector channel due to the 
large exponent ja, the amount of DVs in V + A is basically the same as in V. In 
our opinion, this example shows that one should understand the V and A channels 
individually before drawing any definite conclusion about V + A. 



5 Impact of Duality Violations on R T and a s 

Once the DV functions (13.11) . along with the values for the parameters (I4.2|4.3p . are 
known, one can use Eq. (12.81) to calculate the contribution from our DV ansatz to R T . 
This is given by 

r va = {R v, A) (o) _ 12 nS EW \V ud \ 2 r — 2 (l - -LV (l + ImAvM 

J m 2 m l T \ m l T ) \ miJ 

= «> A )(°) + 5R V T A , (5.1) 

where (R^' A )^ corresponds to the case of no DVs and, consequently, SR^' A corresponds 
to the contribution from duality violations. In the case of our ansatz (13.11) . we obtain 

xdV,a _ k v,a g c , T/ ,2 

SRr ' ~ WUT^J^4 Sew}V - J 

$i(m 2 ) sin (a v ,A + (3v,A m l) + ®2( m l) cos ( a v,A + Pv,A m l) } > ( 5 - 2 ) 
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with 

*i(m?) = 2f3 4 VA - \1l\ A P\ A + Hv,a - (^1v,aPv,a + Ha^va ~ lv, A ) m r > 

$ 2 «) = S^aPvAiIa ~ Pv,a) ~ (Pv,a - H,a0v,a - Hv, A f3v,A)m 2 T , (5.3) 

for s = m 2 . Note that 8R^' A goes to zero for jv,A —> 00 > but is finite in the jv,A — ¥ 
limit. Upon substituting the values from the vector and axial- vector fits (14.2114.31) one 
finds 

8RX = (-0.022 ± 0.013) S EW \V ud \ 2 ; 8R A = (0.0003 ± 0.0003) S EW \V ud \ 2 , (5.4) 

which translates into 8R^/R^ ~ -(1.5 ± 0.9)% and 8R A /R A ~ (0.02 ± 0.02)%. The 
effect of DVs is much smaller in the axial-vector channel, because of the larger 7,4, as 
anticipated. 

For comparison, Ref. [5] takes the resonance model in Ref. [H] as a model for DVs, 
matches it directly to the V + A spectral function near s = m 2 , and finds that the effect 
of DVs is at most 8R^ +A ~ 0.0021,W|Kd| 2 . In our case we keep only the form of the 
asymptotic behavior of this type of model for the spectral function (cf. Eq. (13.11) ). and 
we then fit the parameters of this ansatz to the V and A spectral functions separately. 
We thus find for the effect in V + A the sum of the two results in Eq. (15. 4p . which is 
ten times larger than the estimate of Ref. [5] , but still consistent with the experimental 
data. Ref. [5] also considers an instanton-based model [H]. In this case, the model 
does not have the exponential suppression of Eq. (13.11) . and the contribution becomes 
8R^ +A — Q Ew\Vud\ 2 ■ As it turns out, this is closer to our estimate in magnitude, 
although opposite in sign. 

In order to estimate the shift in a s induced by these values for 8R^' A , we first note 
that the perturbative expression for R^ +A can be written as 

(R V T +A ) PT = N c S EW \V ud \ 2 [l + 8 PT ] , (5.5) 

with 5pt given by the following expansions [5] 

8^ = 0.1082 

0.0334 
8% = 0.1479 

0.0122 

depending on whether one uses the CI or FO prescription!^ Using the experimental 
value for the tau decay ratio R^ +A plus a conservative estimate for the main contri- 
butions from the condensates, Ref. [8] estimates the phenomenological value for the 

5 The last term in these expansions is an estimate of the 0(a b s ) term from Ref. [8]. 



(2^)+ 0.0609 (^£%)V + 

V 0.34 J V °- 34 / 

( 3 + 0.0174 ( 4 + 0.0088 ( «%T) 5 , (5.6) 

V 0.34 J V °- 34 / V 0.34 ) y ' 

(^hl) + 0.0297 + 

V 0.34 J V 0.34 J 

( 3 + 0.0086 ( ^) 4 + 0.0038 ( ^1) ' , (5.7) 

V 0.34 J V 0.34 J V 0.34 J K ' 
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parameter 8 as 

8 phen = 0.2042 ± 0.0050 . (5.8) 

Equating 1 + 8pt + 8R T y + A/ R t ,v+a to 1 + 8 p h en , one obtains an estimate for the shift 
in a s due to the DV contribution (15.21) . yielding 

8a s (m T ) ~ 0.003 - 0.010 , (5.9) 

both for CI and FO prescriptions. The spread of values in Eq. (15. 9p reflects the error 
in the sum of the V and A results in Eq. (15.41) . i.e., in 8R^ +A . 

We consider our result, Eq. (15.91) . as a fair estimate of the systematic error associated 
with duality violations in a s (m T ), as determined from the total nonstrange tau decay 
width. However, the value of a s {m T ) and the condensates may also be determined 
from a combined fit to a set of moments written in terms of pinched weights (see, 
e.g., [H [7]). In this case, our estimate is probably not good enough and a full analysis 
of these moments in the presence of DVs is required, for instance along the lines of 
our recently proposed iterative method described in Ref. [16J. At any rate, it is clear 
that the error associated with DVs is not numerically negligible, at least in the vector 
channel. 

6 Inclusion of e + e~ data 

Since, according to Eq. (12. 9p . DVs entail an extrapolation to high energies and, with 
our ansatz, they turn out to be sizeable in the vector channel, it makes sense to ask 
whether e + e~ data, which extend to higher energies than tau data, can be used to 
further restrict the range of values for the DV parameters (14 .2p . 

However, one cannot relate DVs in tau decay to those in e + e~ without further 
assumptions. Since tau data concerns currents with a flavor structure of the type ud, 
whereas e + e~ data see the flavor-sine let combination^ (2/3)5m-(l/3)dd-(l/3) ss, there 
is an OZI suppressed contribution in e + e~ which is absent in tau decay. Although this 
contribution is 1/N C suppressed and, in perturbation theory, shows up only at 0{a z s ), 
it is not clear how much it might affect the value of the DV parameters appearing 
in Eq. (13.11) . since they encode nonperturbative effects on the time-like axis. As an 
exploratory step, we will simply assume that OZI contributions do not significantly 
affect the values of the parameters in the DV function (13.11) . This is not the only 
difficulty, however. The strange quark mass is much larger than the up and down 
masses, and this should affect the values of the DV parameters for the ss component. 
Fortunately, taking guidance from the model underlying the form of our ansatz [14, 
[T5l \W\ [T7] . one finds that, while the parameters Kv,1v an d fly are related to the 
universal slope of the daughter Regge trajectories, the value of ay is related to the 
mass (squared) of the lowest lying resonance of the Regge tower. This suggests that, 
while the value of ay should be similar in the uu , dd and ud channels (assuming isospin 
symmetry), in the ss channel it should be shifted by a certain amount proportional to 
the strange quark mass. Since we do not know how large this shift may be, in practice 
we introduce a new parameter a' v for the ss component. 

6 We will restrict our discussion to energies below the cc spectrum. 
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Figure 4: Plots of the vector spectral function in tau decay (left panel) and e + e~ (right panel) 
as compared to perturbation theory (green flat line) and the result of the combined fit 16. 2\) 
(red oscillating curve). 

Therefore, in the context of e + e~ data, the spectral function including the ansatz 
for DVs becomes^ 

-ImIT y ( S ) = ^L[l + p( s )] + 

5 1 

- k v e~ lvS sin (a v + (3 v s) + - k v e~ lvS sin (a' v + (3 v s) , (6.1) 

6 6 

where the new parameter a' v takes this aforementioned shift into account, and the 
weights 5/6 and 1/6 correspond to the sum of the squares of nonstrange and strange 
electric charges, after the spectral function is normalized as in tau decays, cf. Eq. ( 14. ip . 

We have attempted a fit of Eq. (16.11) to the e + e~ data[§ in the window 1.1 GeV 2 < 
Smin ^ 9 GeV 2 but, as it turns out, the quality of these data is not good enough to 
allow a stable result for this fit. Another issue of concern is that the e + e~ points come 
from a compilation of different experiments, each with its own systematics, and it is 
not clear to us how to take the systematic errors properly into account (see Ref. [21] 
for a discussion). 

Therefore, we decided to proceed as follows. Since the problem with tau data is that 
they do not extend above s = m 2 , we have supplemented these data with the spectral 
function from e + e~ in the range 4 GeV 2 < s min < 9 GeV 2 as provided by BES [20], at 
present the most precise (inclusive) measurement in this energy range. Inconsistencies 
between exclusive and inclusive determinations of the e + e~ cross section in the region 
3 GeV 2 < s < 4 GeV 2 have prevented us from using these data as well. This is a 
longstanding issue (see, for instance, Ref. [2T]) which, to the best of our knowledge, has 
not yet been fully resolved. Because of this, and because of all the extra assumptions 
we had to make to arrive at Eq. (16.11) . we emphasize that the present exercise should 
be seen only as a very preliminary estimate of constraints from e + e~ data, rather than 

7 Because of OZI-suppressed contributions, the function p(s) in e + e~ differs from that in tau de- 
cay at O(a^), as we have already mentioned. However, the result of the fit is insensitive to these 
perturbatve details, and we will not take them into account. 

8 These data are available at |http: / / pdg. lbl.g ov/ current /xsectTl 
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as a full-fledged analysis. For a more reliable analysis a better understanding of the 
theoretical aspects and better data in the 2 — 4 GeV 2 window would be needed. The 
latter may soon be provided by BaBar [22] . 

With all these caveats in mind, we are now ready to address the question as to 
whether the values of the DV parameters in Eq. (14.21) are consistent with e + e~ data 
above s = m 2 . To try answer this question, we have performed a combined fit in which 
we have taken the ansatz 04. ip to fit the tau data in the region 1.1 GeV 2 < s min < m 2 
and, simultaneously, the ansatz (16.1 j) to fit the e + e~ data in the region 4 GeV 2 < s min < 
9 GeV 2 . The result of this combined fit yields the following values: 

k v = 0.024 ±0.004 

7v = 0.40 ±0.12 

a v = 1.82 ±0.19 

Pv = 2.14 ±0.11 

a' v = 5.2 ±1.4 

X 2 22 



dof 87 ^ 

As one can see, our ansatz produces a reasonable fit to the data with values for the 
parameters which are compatible, within errors, with the results in Eq. (14.21) . There 
is a certain trend towards smaller DVs, in comparison with our results from tau data 
alone. We note that if we compute the x 2 with the parameter values given in Eq. (16.21) 
using only the (nine) BES data points, we find a value x 2 /dof = 10/9. The comparison 
of the fits with the data, and with perturbation theory, can be seen in Fig. HI 

Using the values (16.21) . we can estimate the DVs by substitution into Eq. (15. 2p . 
which yields 

SRV = -(0.0065 ± 0.0042) S EW \V ud \ 2 . (6.3) 

Following the same steps leading to Eq. (15.91) . the corresponding shift in a s , as deter- 
mined from the tau decay width, is given by 

5a s (m r ) ~ 0.001 - 0.003 . (6.4) 

This value is roughly a factor three smaller than that previously found in the fit to tau 
data alone, Eq. (I5.9p . but it is not negligibly small in comparison with the theoretical 
errors quoted in Eq. (II. ip . 



7 Conclusions and outlook 

In this work we have considered the problem of evaluating the possible size of DVs 
in tau decay. Even though there are good reasons to believe that DVs are present in 
QCD, our lack of understanding of these effects makes it very difficult to incorporate 
them in a systematic analysis of the data. However, ignoring the issue altogether may 
lead to a underestimation of systematic effects in the extraction of a s , quark masses 
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and condensates. Some estimate of the effect of DVs is needed, given the high level of 
accuracy presently sought. 

A general framework for the contribution of DVs is provided by Eqs. (12.8112. 91) . In 
order to get to a quantitative estimate of DVs, however, more detailed information 
about the functions ImAy^(s) is needed. In order to do this, we proposed an ansatz 
(13. ip which is based on the models of Refs. [HI [15], [16] [17], but which we expect to be 
of more general validity. Although this validity is speculative, at least it has the virtue 
of being falsifiable with future improvements in experimental data. 

At present, the ansatz is perfectly compatible with the tau data, producing good 
fits for both the vector and the axial-vector spectral functions, and allowing us to 
determine the parameters of the ansatz, cf. Eqs. (14.2114. 3p . These parameter values 
lead to the estimate (15.41) for the contribution from DVs to R T , and, from that, to a 
possible theoretical error in a s (m T ) of order 0.003-0.010. Since we find good fits for 
both vector and axial- vector channels, our fits also give a good description of the V+A 
spectral function. However, we find contributions to R T which are much larger than 
those based on a similar model in Ref . [5] , which were obtained from an estimate based 
on the V + A spectral function only. 

Inclusion of e + e~ data, plus some assumptions, modifies the values of the parame- 
ters in the ansatz to those given in (I6.2p . These modified values lead to a reduction of 
the effect of DVs on R T . This is just a consequence of the flatness of the inclusive e + e~ 
data above 4 GeV 2 [20], to which we have limited ourselves here, as they appear to be 



more reliable than the exclusive data at lower energies. However, given that exclusive 
data below 4 GeV 2 are not consistent with this behavior, it will be important to clarify 
this issue, once more and better data become available. At any rate, even with the 
inclusion of the e + e~ data of Ref. [20J, we find that DVs may lead to a systematic 
error which is as large as 5a s (m T )\ t h — 0.003, i.e., of the order of half the systematic 
errors from all other sources taken together in [5] [8], cf. Eq. (jl.ip . In a conservative 
approach, systematic errors with different origin ought to be added linearly. In light 
of these results, we conclude that future QCD analyses based on tau data should no 
longer neglect the possible contributions from DVs. 
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Appendix: Estimate of OPE contributions to Im Uv,a 



In the analysis of Sec. 4, the analytic form of the ansatz used to fit the spectral function, 
Eq. (14.11) . assumed that the contribution from the OPE to the imaginary part of the 
correlator is negligible. In this appendix we will estimate the imaginary contributions 
coming from 0(a s ) corrections from dimension four and six operators. 
We start by splitting the vector and axial correlators as 



n 



V.A 



tt(0) , n (2) n (4) n (6) 
L1 V,A ' LL V,A ^ LL v A ' 



V.A 



V,A J 



(A.l) 



where the superscripts refer to the dimension of the OPE terms. The perturbative 
contributions are given by p] 



n (0) _ n (o) 

U V ~ LL A 



i . Q 2 



l + 



a s (Q) 



7Y 



(A.2) 



(2) 

Contributions from H VA are proportional to the u, d quark masses and can be safely 
neglected. The contribution from dimension four operators is, from Ref. [1] 
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whereas for dimension six operators one finds [1] 
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(A.4) 



where in all the equations above (• • • ) stands for higher order contributions in the 
strong coupling constant. In the dimension six contribution we have kept only the 
leading log contribution (finite pieces will not contribute to the imaginary part), and 
we used the following short-hand notation for the four-quark condensates: 



{ij>^T a W'fT a 1>) , 



Assuming vacuum factorization [2], the four-quark condensates simplify to 



(A.5) 



2 , T 



In order to calculate ImHy^t), one uses that 



a s (Q) = a s {/j) <; 1 + — — ^— log -j- 



with 



As a result, 



2ra/ - 11JV C 



6 



rtf=3 



9 
2 



Im[a s (Q)] = i [« s (g)] 2 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



It is now straightforward to calculate the different OPE contributions to the spectral 
functions. For normalization purposes it is convenient to work with v\ = 27rlm Hy, a\ = 
27rIrnII J 4 . Taking a s (l.l GeV 2 )/7r ~ 0.15, one finds for the perturbative contribution: 



(0) 



oc s {Q) 



71 



0.58 . 



(A.10) 



Taking {a s G 2 } ~ 0.017T GeV 4 , the gluon condensate contribution amounts to 



Svf 2) = 5a^> ~ -4 • 10 



(G 2 ) 



GeV 



(A.n; 
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For the quark condensate contribution, we will assume SU(3) symmetry in the conden- 
sates, the PCAC relation (m u + mj(^) = —m%f%, and estimate m s /(m u + m^) ~ 10. 
This yields 

5v m» „ _ 3 . 10 -4 ^ } Sa m)) „ 1 . 10 -4 (A12) 

In Ref. pQ it was argued that the biggest non-perturbative contributions to R T come 
from the dimension six operators. There is some evidence that naive factorization does 
not work for these condensates. Based on a phenomenological fit, Ref. [1] concluded 
that there might be an enhancement over the result from factorization that can be 
represented by rescaling the condensate 

<^) 2 — A(^> 2 , (A.13) 

where A is a fudge factor which, for typical values of a s and the quark condensate, can 
be estimated to be A ~ 3 — 9. We take the quark condensate (^V) = (— 240MeV) 3 
and, to be conservative, we take A = 10. Accordingly, 

^ M) <4-10-^, ia ?W <_ 10 -* . (A . 14) 

Despite the enhancement by one order of magnitude, the contribution of dimension- 
six operators (as well as the dimension-four operators) is negligible compared to the 
perturbative contribution, Eq. (1A.10I) . We conclude that we can safely neglect the con- 
tribution of the OPE to the spectral functions, in comparison to perturbation theory. 
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